Rapid Spin Deceleration of Magnetized Proto-Neutron Stars 
via Asymmetric Neutrino Emission 



Tomoyuki Maruyama, 1 ' 2 Jun Hidaka, 2 Toshitaka Kajino, 2,3 Nobutoshi Yasutake, 4 
Takami Kuroda, 2 Myung-Ki Cheoun, 2,5 Chung- Yeol Ryu, 6 and Grant J. Mathews 7 
College of Bioresource Sciences, Nihon University, Fujisawa 252-8510, Japan 



o ■ 

£N| ■ 2 National Astronomical Observatory of Japan, 

C 1 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan 

[ 3 Department of Astronomy, Graduate School of Science, 

. University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo 113-0033, Japan 

* Department of Physics, Chiba Institute of Technology, 

2-1-1 Shibazono, Narashino, Chiba 275-0023, Japan 

_J , Department of Physics, Soongsil University, Seoul, 156-743, Korea 

Q-H, 6 General Education Curriculum Center, Hanyang University, Seoul, 133-791, Korea 

7 Center of Astrophysics, Department of Physics, 
& ' ' rn» 

C/3 ' University of Notre Dame, Notre Dame, IN 46556, USA 

a ' 

1 'I (Dated: February 1, 2013) 

, We calculate the spin deceleration of proto-neutron stars by asymmetric neutrino absorp- 

m ■ 

■ tion in the context of a fully relativistic mean field theory. We calculate for the first time the 

' spin deceleration of neutron stars due to asymmetric neutrino absorption in a toroidal mag- 

| netic field configurations. We find a surprising effect that the deceleration can much larger 

, for asymmetric neutrino absorption in a toroidal magnetic field than the usually presumed 
braking due to magnetic dipole radiation. This may explain the observation that magnetars 

> ■ 

. ' appear to have had a more rapid deceleration in the past. 

J] , PACS numbers: 25.30.Pt, 26.60.-c, 24.10. Jv, 97.60.Jd 



Magnetic fields play an important role in many astrophysical phenomena. The observed asym- 
metry in supernova (SN) remnants, pulsar kick velocities and the existence of magnetars j^Q] 
all suggest that strong magnetic fields affect the dynamics of core-collapse SN explosions and 
the velocity 4j that the proto-neutron stars (PNS) receive at birth. 

There are now two major supernova explosion scenarios leading to asymmetric morphologies 
in observed supernova remnants. One is the standing accretion shock instability (SASI)-aided 
neutrino driven explosion [|], Q|. The other is the magneto-rotational explosion (MRE) s|- 
Both mechanisms may also be a source for pulsar kick velocities |{], [h]]. The MRE takes place 
through the extraction of the rotational energy of the PNS via strongly amplified magnetic 
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fields ~ 10 15 G [39]. Thus, it is expected to leave a magnetar remnant. However, there are many 
unknown aspects of these scenarios, such as the progenitor and final core rotation and magnetic 
field profile. Hence, there is not yet a definitive understanding of the observed asymmetry and 
remnant kick velocities in core collapse supernovae. 

Moreover, there is accumulating evidence that magnetars spin down faster than expected [jjj. 
Usually it is supposed that magnetars spin down mainly by dissipating rotational energy into 
magnetic dipole radiation. However, it is possible that they could have a much more rapid spin 
deceleration early in their formation. A new source for this deceleration is the main point of the 
work reported here. 

Here, we show that the early spin deceleration of a PNS could result from an asymmetry in 
the neutrino emission that arises from parity violation in weak interactions 
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an alternative scenario see Ref. 



14( and/or an 



asymmetric distribution of the magnetic field [151 ] in strongly magnetized PNSs. (However, for 



161].) Recent theoretical calculations 







, [li| have suggested 



that as little ~1% asymmetry in the neutrino emission out of a total neutrino luminosity of 
~ 10 53 ergs is enough to explain the observed pulsar kick velocities. Here we show that such 
emission could also explain a rapid early spin down rate of magnetars. 

Both static and dynamic properties of neutron-star matter have been studied 19h21|] at high 



temperature and density in the context of spherical non-magnetic neutron star models. Such as- 
pects as an exotic phase of strangeness condensation 
powered thermal evolution 
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241 ] . nucleon superfluidity [251 ] . rotation- 



, a quark-hadron phase transition 



27] ) etc., have been considered. 



Neutrino propagation has also been studied for PNS matter including hyperons 



28(]. These the- 



oretical treatments of high-density hadronic matter, however, have not yet considered the effects 



of strong magnetic fields. Although previous work 1t], [l8( has studied the effects of magnetic 
fields on the asymmetry of neutrino emission, the neutrino-nucleon scattering processes were 
calculated in a non-relativistic framework 



n 



171 ] and only a uniform dipole field configuration was 



considered. 

In our previous papers 
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301 ] we discussed the neutrino scattering and absorption cross sections 
in hot and dense magnetized neutron-star matter (including hyperons) in fully relativistic mean 
field (RMF) theory 32j. Those papers demonstrated that magnetic fields enhance the scattering 
cross-sections in the direction parallel to the magnetic field for neutrinos while also reducing the 
absorption cross-sections in the same direction. When the direction is anti-parallel, the opposite 



[39] A few processes can be candidates for the amplification mechanisms such as winding effect or the magneto- 
rotational instability fl"ll |. 



occurs. 
It was shown in Ref. | 

B ~ 2 x 10 17 G, the enhancement of the scattering cross-sections is ~ 1% at a baryon density 
of pb = 3/3o, while the reduction in the absorption cross section is ~ 2%. This enhancement 
and reduction were shown to increase the neutrino momentum flux emitted along the direction 
of the magnetic field and to decrease the emitted momentum flux emitted antiparallel to the 
magnetic field. This asymmetry was then applied to a calculation of pulsar-kick velocities in 
the context of a one-dimensional Boltzmann equation including only the dominant effect of 
neutrino absorption. PNS kick velocities of ~ 610 km-s -1 were estimated for a neutron star of 
baryon mass M^s = 1.68Mq, a uniform dipole magnetic field of B = 2 x 10 17 G, an isothermal 
temperature of T = 20MeV, and a total emitted neutrino energy of Ex ~ 3 x 10 53 erg. 
However, magneto-hydrodynamic (MHD) PNS simulations (e.g. [jy]]) have shown that the mag- 
netic field inside a neutron star may have a toroidal configuration . Here, we apply the above 
results on asymmetric neutrino absorption to this case. One can easily suppose that neutrinos 
are preferentially emitted along a direction opposite to that of the rotation. This would enhance 
the spin down rate of a PNS. In this letter, we present for the first time a study of the effect of 
asymmetric neutrino absorption on the spin deceleration of PNSs. 

Even a strong magnetic field has less mass-energy than the baryonic chemical potential in 
degenerate neutron-star matter, i.e. \J eB <C p^, where p\, is the baryon chemical potential. 
Hence, we can treat the magnetic field as a perturbation. We then ignore the contribu- 
tion from convection currents and consider only the spin-interaction. We also assume that 



\pbB\ <C E£(p) = yp + M£ , and treat the single particle energies and the wave function in a 
perturbative way. 

In this framework we then obtain the wave function in a magnetic field by solving the Dirac 
equation: 

- Ml - U (b) l0 - p b Bo z ] u b (p, s) = 0, (1) 

where M£ = M b — U s (b), and U s (b) and Uo(b) are the scalar mean- field and the time-component 
of the vector mean-field for the baryons b. These scalar and vector fields are calculated in the 
context of RMF theory. 

In Rcfs. we calculated the absorption cross-section a a in PNS matter for an interior 

magnetic field strength of B = 2 x 10 17 G and a temperature of T = 20 MeV. In those results 
we demonstrated that the absorption cross-sections are suppressed in the direction parallel to 
the magnetic field B by about 2 — 4% in the density region of pB = (1 — 3)/3q- The opposite 
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effect occurs in the anti-parallel direction. The net effect of these changes in the absorption 
cross sections leads to an increase in the emitted neutrino momentum flux along the direction 
of the magnetic field and a decrease of the momentum flux emitted in the antiparallel direction. 
However, it is quite likely 



311 ] that the magnetic 



field exhibits a toroidal configuration within 
the PNS. Hence, we now consider the impli- 
cations of a toroidal field configuration on neu- 
trino transport in a strongly magnetized PNS. 
For this purpose we solve for the neutrino 
phase-space distribution function f u (r,k) us- 
ing a Boltzmann equation as described below 



and in Rcf. 



30]. 
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We assume that the system is static and nearly 
in local thermodynamic equilibrium. Under 
this assumption the phase-space distribution 
function satisfies df u /dt = and can be ex- 
panded as f u (r, k) = fo(r, k)+Af(r, k), where 

c . , , . n FIG. 1: (Color online) The upper panel (a) shows 

the first term is the local equilibrium part, and 




the second term is its deviation. 
Furthermore, we assume that only the domi- 



the baryon density distribution for a PNS with T = 

20 MeV and Y L = 0.4. The solid and long-dashed 

lines show results with and without As, respectively, 
nant effect of absorption contributes to the neu- ThfJ mkMle pand (b) showg ^ numbcr fractions 

trino transport, and that the neutrinos travel f or pro tons, lambdas and neutrinos in a PNS. The 
along a straight line. The Boltzmann equation solid and dashed lines show the proton fraction in 



for f v then becomes 
- d 



dr de 
a A {r,k 



V 



dr 

A/(r,fc), (2) 



systems with and without As. The dot-dashed line 
indicates the lambda fraction. The dashed and dot- 
ted lines denote the neutrino fraction in systems 
with and without As, respectively. The lower panel 
(c) shows the field strength distribution at z = for 



where e„(r) is the neutrino chemical potential the toroidal magnetic fields considered here. The 



at coordinate r. 

Here, we define the variables xl = (r ■ k)/\k\ 
and Rt = r — (r ■ k)k/k 2 . We can then solve 
Eq. ([2|) analytically 



solid and dashed lines represent those for ro = 8km 
(Mag- A) and 5km (Mag-B), respectively. 



Af(x L ,R T ,k) 



■•-L 



dy 



de u df 
dx L de u 



cxp 



B * a A (z,R T ,k) 
az y 



(3) 
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where the center of the neutron-star is at r = (0, 0, 0), and all of the integrations are performed 
along a straight line. 

In this work we utilize an equation of state (EOS) at a fixed temperature and lepton fraction 

n nn 

by using the parameter set PM1-L1 [33[] for the RMF as in previous work |29|, 1301 ] . 

We show the baryon density in Fig. [T^t and the particle fractions in Fig. [lb as a function of the 

neutron-star radius. For this figure we assume a neutron star baryonic mass of Mns = 1.68Mq, 

a temperature of T = 20MeV, and a lepton fraction of Yl = 0.4. The moment of inertia of the 

neutron star becomes Ins = 1-54 x 10 45 g-cm 2 or 1.36 x 10 45 g-cm 2 with or without lambdas, 

respectively. 

One can see in Figs. [T^i and [lb the appearance of lambdas for a baryon density greater than 
about twice the saturation density of nuclear matter, i.e. ps > 2po> where po « 2.7 x 10 14 g 
cm -3 . This softens the EOS and leads to an increase in the baryon density and neutrino fraction 
for r < 8km relative to hadronic matter without lambdas. 

The ratio of the total rate of angular momentum loss to the total power radiated by neutrinos 
at a given spherical surface Sn is 

(cdL z /dt\ _ Is N dn rf^^f {r,k){r x k) ■ n 

fs N f (f^ 

where n is the unit vector normal to Sn- For illustration we will consider surfaces for which 
p B = Po, and p B = po/10. 

In this model the neutrino current is assumed to be constant. Hence, we can use the chain rule 
to write dA/dt = cdA/dxi {A = L z ,Et) on the boundary surface, Sn- We can then obtain the 
angular acceleration from the neutrino luminosity, C u = (dEr/dt), 

d NS \ dEr/dt) 

For a NS spin period P, the angular velocity is u = 2tt/P, and the angular acceleration is 
-2irP/P 2 . Thus, we obtain 

P _ P ( cdL z /dt \ 

P 2ircI NS \ dEr/dt J v ' U 
We adopt the following parameterization of the toroidal magnetic field configuration, 

B = B G L (z)G T (r T )e^ (7) 

where = (— sin <f>, cos 4>, 0) in terms of the azimuthal angle <f>, and 

G L (z) = — G T (r T ) = — 7 . (8) 

M + e z l ao ] [l + e( r 'T- r o)/«ol z 



\ dEr/dt J J s dn r J^Af(r,k)k.n ' 
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For our purposes we assume that the toroidal magnetic field is aligned along the direction of 
the spin rotation. Toroidal fields could however be oppositely oriented and can even invert with 
time. The generalization to other orientations would imply that neutrino emission could also 
accelerate the spin of the star. Another complication is the realistic case of a double torus 
configuration for the magnetic fields. Assuming that the directions of the toroidal magnetic 
field and the spin rotation are the same in the northern hemisphere, our effects would brake and 
decelerate the rotation. In the southern hemisphere, however, the direction could be anti-parallel 
and the asymmetry in neutrino absorption may even accelerate the rotation. This might lead 
to a complicated twisting mode that could induce p-mode oscillations. 

In Fig. [lb we show the magnetic field strength \B/B$\ for different field configurations, with 
ao = 0.5 km and tq = 8.0 km (Mag-A) or tq = 5.0 km (Mag-B). We assume typical values for 
the neutrino luminosity C u ~ 3 x 10 52 erg-s -1 [l8] and a magnetar spin period of P = 10 s 
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35f | . We use B = 2 x 10 17 G and summarize our results in Table HI 





Comp. 


dcL z j dt 


P/P (s- 1 ) 


dE T /dt 


PB = Po PB= Po/10 MDR 


Mag-A 


p, n 


66.9 


7.82 x 1CT 2 1.03 x 10~ 2 7.76 x 10~ 12 


p, n, A 


109 


1.13 x lO" 1 1.11 x 10~ 2 9.86 x 10~ 12 


Mag-B 


p, n 


9.64 


1.13 x 10~ 3 4.50 x 10~ 5 7.76 x 10~ 12 


p, n, A 


7.81 


8.07 x 10~ 4 2.29 x 10~ 5 9.86 x 10~ 12 



TABLE I: The 1st column gives the model for the toroidal magnetic field configuration (see text). The 
2nd column shows the presumed composition of nuclear matter, i.e. "p, n" for nucleonic and "p, n, A" 
for hyperonic matter. The 3rd column denotes results from Eq. the 4th and 5th columns are results 
obtained using Eq. ([5]) at the indicated baryon density. The 6th column shows the usually presumed 
spin-down rate from magnetic dipole radiation [Eq. ([9])]. 



To illustrate the density dependence of this effect we include the results obtained for Sn located 
at pb = po/10 in the fifth column. 

For reference the last column of Table [J shows P/P calculated with the usual magnetic dipole 
radiation (MDR) formula 361 ] . 



PP = B' 



( 3I NS c 



(9) 



y 1 2 ( ui^ y 1 

V 8tt 2 ^6 J \\2^ 2 I 2 NS ) 

Here R is the NS radius whose value is related to the moment of inertia, Ins = 2MnsR 2 /§■ 
Since MDR derives from the external magnetic field we utilize a typical exterior magnetar 
poloidal magnetic field of B = 10 15 G in Eq. ©. From this we see that the spin deceleration 
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from asymmetric neutrino emission can be much larger than that for MDR when the neutrino 
luminosity is high. 

In this calculation we have ignored the neutrino scattering and production processes. The 
neutrino scattering process enhances the asymmetry of the emission, although its contribution 
to the mean-free path is much smaller than that from absorption in the density region of interest, 



Po<PB< 3p 



Neutrino production in a magnetic field is known to cause asymmetry in the neutrino emission 
. The cross-section for the neutrino production reaction, e~+p — > n, A+f e , is qualitatively 
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the same as that for the absorption reaction, v e + n, A — > p + e~ . The only difference is the 
small contribution from the magnetic part of the initial and final electron states. Hence, this 
production process would tend to enhance the asymmetry and to also contribute to the spin 
deceleration. 

Even so, this estimation of P/P still has some uncertainties due to unknown variables such as 
the interior strength of the magnetic field, the spin period of the NS core, etc. Nevertheless, 
since our value of P/P is at least 10 5 times larger than that for MDR, the asymmetric emission 
could be significant under many possible early conditions. Moreover, other processes such as 
the neutrino scattering and production tend to increase the spin deceleration. Thus, we can 
conclude that asymmetric neutrino emission from PNSs may play a very important role in the 
spin deceleration of a magnetic PNS. 

In summary, we have estimated the spin down of a PNS with a toroidal magnetic field con- 
figuration by considering the asymmetric neutrino absorption cross-sections calculated in the 
context of RMF theory. We then solved the Boltzmann equation using a one-dimensional at- 
tenuation method, assuming that the neutrinos propagate along an approximately straight line, 
and that the system is in quasi-equilibrium. We only included neutrino absorption which dom- 
inates [^j over scattering in producing asymmetric momentum transfer to the PNS. We found 
that asymmetric neutrino emission can have a very large effect on the early spin deceleration 
of a PNS. Indeed, this effect can be much larger than the braking from a magnetic dipole field 
configuration. 

From Eq. ([6]) we can see that P/P from asymmetric neutrino emission is proportional to P, 
while P/P oc P~ 2 for MDR. This fact implies that asymmetric neutrino emission could play a 
significant role well into the later stages of core-collapse supernovae. 

Finally, we caution that definitive conclusions should involve a fully dynamical MHD simulation 
of the evolution of a PNS with asymmetric neutrino scattering and production as well as ab- 
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sorption in a strong magnetic field. Nevertheless, the results presented here make the case for a 
significant influence of asymmetric neutrino absorption on the process of magnetar formation. 
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